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ABSTRACT 

We present a detailed analysis of the Chandra High Energy Transmission Grating Spectrom- 
eter (HETGS) and XMM-Newton high resolution spectra of the bright Seyfert 1 galaxy, Mrk 
290. 

The Chandra HETGS spectra reveal complex absorption features that can be best de- 
scribed by a combination of three ionized absorbers. The outflow velocities of these warm 
absorbers are about 450kms~ 1 , consistent with the three absorption components found in 
a previous far UV study. The ionizing continuum of Mrk 290 fluctuated by a factor of 1.4 
during Chandra observations on a time scale of 17 days. Using the response in opacity of 
the three absorbers to this fluctuation, we put a lower limit on the distance from the ionizing 
source of 0.9 pc for the medium ionized absorber and an upper limit on distance of 2.5 pc for 
the lowest ionized absorber. The three ionization components lie on the stable branch of the 
thermal equilibrium curve, indicating roughly the same gas pressure. Therefore the thermal 
wind from the torus is most likely the origin of warm absorbing gas in Mrk 290. 

During the XMM-Newton observation, the ionizing luminosity was 50% lower compared 
to the one in the Chandra observation. The RGS spectrum is well fitted by a two phase warm 
absorber, with several additional absorption lines attributed to a Galactic high velocity cloud, 
Complex C. Neither the ionization parameter £ nor the column density Nh of the two absorb- 
ing components varied significantly, compared to the results from Chandra observations. The 
outflow velocities of both components were 1260 km s -1 . We suggest that an entirely new 
warm absorber from the torus passed through our line of sight. 

Assuming the torus wind model, the estimated mass outflow rate is ~ IMq per year, 
while the nuclear accretion rate is <~ 0.5M Q per year. The O VII and Ne IX forbidden lines 
are the most prominent soft X-ray emission lines, with a mean redshift of 700 km s" 1 relative 
to the systematic velocity. There seems to be no relation between emission lines and warm 
absorbers. 

Key words: galaxies: Seyfert — galaxies: absorption lines — X-rays: galaxies — galaxies: 
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1 INTRODUCTION 



Ionized outflowing gas (or 'warm absorber', WA) is believed to be common in the in ner parts of Active Galactic Nuclei (AGN) based 
on X - ray spectroscopy. About 50% of Seyfert Is i Reynolds! 1997 ; George et al. 1998 ) and quasars l Piconcelli et al. 2005 



Misawa et al 



2007 ; Ganguly & Brothertonl 2008 ) are characterized by blue-shifted X-ray absorption lines, the velocities o f which are usually se veral 
hundred kms" 1 . UV narrow absorption lines from outflowing gas were also detected at a si milar fraction i Crenshaw et al. I ll999h . and 



showed velocities and ionization states as similar to WAs in some cases dGabel et al.ll2005al lbl: I Costantini et al 
cases there is no clear physical correspondence between them i KrissfcoO^ ). 



J I2OO7I) 



. However, in other 



Several stu dies have shown that WAs can be well modeled by including mul tiple absorbing components, which appear to be in pressure 

Krongold et al. 200S ). Thus , the structure of WAs can be interpreted as a multi- 



equilibrium (e.g. Netzer et alj|2003 



McKernan et al. 2007 



phase medium, and most li kely occupies the s ame spatial region iKrolik & Krissll200lh . Ano ther suggestion that WAs oc cur in continuous 
density-stratified streamers 1 Behar et al]|2003 ; Steenbrugge et alj|2005 ) seems to be ruled out 1 Krongold et al.12005 , 2007 ). 

Unfortunately, very little is still known about the physical states of WAs, or their dynamics and geometry. The most fundamental 
question is where WAs originate. P roposed locations span a wide range in radial distanc e from the central ionizing sou rce: the accretion disk 
( Elvislboodl ICron gold et alj|2007h , the broad line region ( BLR) dKraemer et alj|2005h , he putative obscuring torus dKrolik & KrissluOOll ; 
Blustin et al.l2005 |) and the narrow line region 1 Behar et al|2003; Crenshaw et al.120091) . Simulatio ns have offered similar p ossibilities, such 



as the accretion disk wind 



Proga & Kallman 200 



(P 



34; 



Risaliti & Elvis 



2009 ), the wind from torus 1 Dorodnitsvn et al. 2008J), the large-scale 



outflows jKurosawa & Progall2009l) . but ruled out few options. 

In principle, these outflows could potentially pr ovide a common form of AGN feedback. But whether the o utflow winds cou ld extend 
into the interstellar medium and intergalactic medium ( Nicastro et alj2008 ) or just fall back into the accretion disk 1 Stoll et al .120091) , depends 
critically on the exact location of WAs when estimating the mass outflow rate and the to tal outflow kinetic en ergy. One key diagnostic for the 
location of WAs is based on applying non-equilibrium models on time variable AGNs 1 Nicastro et al.ll 1999h . Such strong variability needed 
to apply to this method is not that common. 

Mrk 290, a relatively nearby Seyfert 1 galaxy (z = 0.0304 ± 0.001 4 from the Sloan Digital Sky Survey), has been extensively studie d 
by low resolution instruments (e.g. Wood et aDll984 ; Kruper et al.lll99of) . The WA was found in the ASCA observation 1 Turner et al.lll996h . 
The source flux of 2-10 keV band was 8.4 x 10" 11 ergs" 1 cm" 2 in the ASCA observation. In 2003, the first high resolution observation 
was carried out using Chandra HETGS, simultaneously with Far Ultraviolet Spectroscopic Explorer (FUSE ) to determine the ionization 
state and kinematics of the WA, and assess association with the UV-absorbing gas. Mrk 290 was observed by XMM-Newton in 2006. 

We present the analysis of the ionized absorbers in Mrk 290, using the high resolution X-ray spectra. This paper is organized as follows. 
In section 2, we describe the data reduction. Section 3 is devoted to the analysis of the data, and section 4 to the detailed model of the ionized 
absorbers of Mrk 290. In section 5 we discuss our results and present our con clusions in section 6. The cosmological parameters used are: 



Ho = 73 km s 1 Mpc , Q.a = 0.72 and Q m = 0.24 l lSpergel et al .1120071) . The luminosity distance is 123.7 Mpc in this cosmological 



model. The quoted errors refer to 90% confidence for one interesting parameter. 



2 OBSERVATION AND DATA REDUCTION 
2.1 Chandra Data: 2 Flux States 



Chandra performed four observ ations of Mrk 290 for a total duration of 251 ks, using the High Energy Transmission Grating Spectrometer 
(HETGS Jcanizares et al]2005h . Table[j]lists the observation log. 



The HETGS is composing of two grating sets, the High Energy Gratings (HEG) and the Medium Energy Gratings (MEG), the resolutions 
of which are 0.012 A and 0.023 A, respectively. All observations were reduced uniformly and in a standard way using the Chandra Interactive 
Analysis of Observations (CIAO) software (Version 4. 1.1). Flux calibration of each observation was carried out using the Chandra Calibration 
Database (Version 4.1.2). The +lst and -1st orders of each HEG and MEG spectra were combined; higher orders were ignored. 

We generated the light curves in 1000 second bins for the short wavelength band (2-10 A), as shown in Fig.Q] The comparison between 
them indicates that the time-average count rate increased by a factor of ~ 2, between observations 4399 and 3567, only half a m onth apart. 
Smaller variations are observed on time scales of 10 4 s within each observation, which were also found bv lTurner et alj 1 1996t) . Figure Q] 



shows the softness ratio, defined as the count rate in the 15 -25 A band divided by that in the 2-10 A band, to illustrate the different variability 
characteristics below and above 1 keV 1 Netzer et alj|2003 ). The softness ratio increases from 0.075 to 0.100 along with a higher count rate, 



between observations 4399 and 3567. In light of these results, we have divided the entire dataset into two groups with higher (ID 3567, 4441 
& 4442) and lower (ID 4399) count rates. 



2.2 XMM-Newton Data Reduction 

Mrk 290 was observed four times by XMM-Newton 1 Jansen et alJbOQl ) during the period April 30 to May 6, 2006. The European Photon 
Imaging Camera (EPIC) was operated in small window mode. Both RGS (Reflection Grating Spectrometer) cameras were in spectroscopic 
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mode. Raw data from all observations were processed with the standard Science Analysis System (SAS) using version 8.0 and the most 
updated calibration files. 

High 'flaring' background periods were all filtered to maximize the signal-to-noise ratio (S/N). While the total exposure time is 100 ks, 
after this screening process, the net exposure time was 65.5 ks, as shown in Table [2] The four XMM-Newton data were merged, since no 
significant flux variability was found. 

The pn spectra were extracted using pattern 0-4. The 0.5-10 keV count rate, 2.5 counts s _1 , was under the pile-up threshold, verified 
with the EPATPLOT. Each source spectrum was extracted from a circular region of 40" of radius centered on the peak of the X-ray emission. 
The background spectrum was taken from a vicinity source-free circular region of radius of 55". RGS data were processed with the task 
RGSPROC according to the standard method proposed in the SAS thread. 



3 SPECTRAL ANALYSIS 



We used the In teractive Sp ectral Interpretation System (ISIS version 1.5.0, dHouckll20020 ) to fit the spectra. When fitting the data, we used 
the C statistic l Cashll 19760 to find the best-fitting model parameters in the high-resolution spectra. 



3.1 The Broadband Continuum Modeling 

The Mrk 290 spectrum is affected by a relatively low Galactic absorption: Nh = 1-76 x 10 cm - l IKalberla et al .120050 . which is included 
in all models. 

The Chandra spectra were binned to 0.01 A, and fitted in t he wavelength range 1.5-25 A, as shown in Fig. [2] In cont rast to previous 
results, where Mrk 290 was described by an absorbed power law i Turner et alJl996tlshinozaki et alfeoodJjiang et alj|2006Q , we found that 
a black-body component was needed in addition to the absorbed power law in order to fit the continuum. The black-body component was 
included to account for the soft excess, that could come from a standard optically-thick accreti on disk. The intrinsic cold absorpti on was 
negligible compared to the Galactic column density. In fact, half of AGN dMcKernan et al.ll2007h and quasars Jpiconcelli et aljEo05l) do not 
appear to have significant intrinsic cold gas absorption. 

The pn spectrum was grouped to give a minimum of 50 counts per bin in order to apply x 2 statistics. The same model as fitted to HETGS 
spectra but with free parameters provided a good fit to the pn data, except for two prominent residuals around 0.8 and 6.3 keV. These two 
features are consistent with typical absorbing features of WA and Fe Ka emission line, respectively. Again, no cold absorption component 
above Galactic was needed. Fig.[3]shows the spectrum in the 0.5-10 keV band. 

The best fit parameters of all spectra are given in Table [3] including the fluxes and luminosities in the 2-10 keV band. The hard X-ray 
luminosity increased by a factor of 1.4 in 17 days between Chandra observations. When XMM-JVewton re-observed the source three years 
later, the luminosity showed a decrease of 50% compared to the highest value measured in the Chandra spectra. So hereafter, we refer to 
these spectra as 'high-', 'mid-' and 'low-state' spectra as labeled in Table[3] 



3.2 Intrinsic Absorption Line Properties 

The HETGS spectra, which are binned to 0.01 A, reveal many absorption line features, some of which are blends of several lines. H-like 
and He-like lines of O, Ne, Mg, Si, S, and L-shell lines of Si, Fe are marked in Fig. [2] The M-shell complex of iron, called the unresolved 
transition array (UTA) is observed as a depression at 16-17.5 A, although it is very shallow in the present data. We identified one prominent 
narrow absorption line of Fe X in the UTA (see Table|4j. All statistically significant lines with Poisson probability greater than 0.995 were 
measured by fitting a negative Gaussian to the continua within ±0.5 A band, as shown in Table|4] The widths of the Gaussian lines were all 
fixed to 1 mA. The outflow velocities about 450 km s~ were derived with respect to the systemic redshift 0.0304, except those lines labeled 
as Galactic absorption.There are fewer significant lines in the mid-state spectrum than in the high-state spectrum. 

The RGS spectra are binned to 0. 1 A, in which the iron UTA shows clearly. We fitted RGS 1 and 2 simultaneously in the wavelength 
band 7-35 A, but for clarity added them in Fig.|4] The continuum model was fixed to that obtained from pn spectrum, except for allowing the 
black body component to vary. We found several highly ionized features: Ne IX combined with Fe XIX, Ne X Lya, C VI, N VII, Fe XVIII and 
Fe XX complex. A Gaussian was fitted to each feature and the results were listed in Table|4] The outflow velocities of these lines are more 
than 1000 kms -1 , differing by the several hundred kms -1 from those derived from the HETGS spectra. 



3.3 Galactic Absorption 

The spectra of Mrk 290 show narrow absorption lines at zero redshift in Table [4] O I absorption at 23.5 A in Fig. |4] is attributed to the 
ISM of our Galaxy ' Wilms et al]|20000 . We also detected ionized material traced by O VIII in the HETGS spectra, and O VII, N VII, N VI, 
and possibly C V in the RGS spectrum. This ionized absorber may be associated with the high-velocity cloud complex C, w hich is a low- 



metal licity gas cloud pl unging toward the d isk and beginning to interact with the ambient gas that surrounds the Milky Way l lCollins et al 



20070 . In the UV band. lPenton et al.1 hoodt) detected absorption lines from complex C using HST observations. Since this ionized gas does 
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not affect the continuum shape l lCostantini et alj|2007l) . we fitted these lines using Gaussians. The upper limits of the infall velocity with 
respect to our Mi lky Way were meas ured from these lines less than fOOO km s _1 , consistent with the values from the UV lines: 200 kms -1 
and 120 kms" 1 JCollins et al]|2007h. 



3.4 Emission Lines 

In the soft X-ray band, Mrk 290 exhibits a prominent Ne IX forbidden line in the high-state spectrum and a O VII forbidden line in the 
high- and low-state spectra. They are red-shifted at offset velocities of 700 kms" 1 with respect to the systemic redshift, which is definitely 
different from the absorption lines. The velocities do not change over three years suggests that these emission lines may be emitted in the 
same gas. The fit results for the three lines are shown in Table [5] 

The iron emission line at 6.4 keV (or 1.9378 A) is present in both the high-state HETGS spectrum and the low-state pn spectrum. This 
fluorescent line is Ka transitions from neutral iron, consistent with Fe I-XII. A Gaussian line, the width of which was fixed at 1 mA, had 
been included to fit the Fe Ka line in the 5-8 keV band. Table[5]lists the measurements. The equi valent widths (EW s) of 48 ± 17 eV for the 
high-state spectrum and of 23 ± 17 eV for the l ow-state spectrum are c onsistent w ith the result of | Shu et al. However they are much 



smaller than the EWs from ASCA data: 273 eV Jshinozaki et alj2006h or 500 eV hurner et aljl99&h 



We allowed the line width to float and fitted the iron Ka again. The Cash statistics changed by 2 and 8 when reducing one d.o.f. for the 
high- and low-state spectra, respectively, and the equivalent widths (EWs) increased to 71 ± 35 eV and 81 ± 29 eV. The FWHM of Fe Ka 
in the high-state sp ectrum is 57OOI4500 kms" 1 , and in the low-state pn spectrum is MOOOl^ooo kms" 1 , which agrees with ASCA result 
1 Turner et al.lll99&h . 



In the high-state spectrum, a line at 1.897 A seems consistent with the n=2-l reson ance transition of Fe XXV (at a rest energy of 6.696 



keV). This line could be produced in Compton-thin, low velocity photoionized material jBianchi et alj|2005 ) 



There are unidentified emission features with Poisson probabilities less than 0.001 in the high-state spectrum at 2.666 A and 8.725 A, 
of which the fluxs are about few x 10 6 phcm -2 s _1 . However, this two features are probably spurious. 



4 PHOTOIONIZATION MODELING 

The publicly available photoionization code XSTAF0 is the tool we used to model the physical conditions of the absorbing gas. It includes 
H- and He-like transitio ns of most ions, L- shell transitions of Fe, and important M-shell transitions of Fe that contribute to the UTA feature. 
We used PVM_XSTAR jNoble et al.ll2009h to g enerate grids of mode ls, which provides parallel execution of XSTAR 2. llnl 1 . 

We adopted the default solar abundances 1 Grevesse et al. 1996b in XSTAR, and used a velocity turbule nce (fc-value) of 170 km s -1 , 
which corresponds approximately to the limiting MEG spectral resolution (~ 300 kms" 1 FWHM at 0.5 keV) dMcKernan et al }200i) . This 
choice is justified by the absence of detectable line broadening. 

The photoionization model used here contains three parameters: red-shift z; total neutral hydrogen column density Nh', and the ioniza- 
tion parameter £ = Li 0n /(n e R 2 ), where Li an is the ionizing luminosity in the range 1-1000 Ryd, n e is the electron density, and R is the 
distance of the ionized gas from the central ionizing source. 



4.1 The Spectral Energy Distribution 

First of all, we constructed the spectral energy distributions (SEDs) of Mrk 290 to generate the photoionization models. However, the intrinsic 
X-ray spectrum, an important part of the SED, is obtained by fi tting the X-ray data w ith photoionization models. Thus a two-step iterative 



process is employed to set the intrinsic spectrum, as described in lYaqoob et al. U2003h 



In the first step, the ionizing X-ray spectrum is obtained by fitting the continuum without a WA to construct a SED. The SED is then 
used to generate a grid of warm absorption models that are fit to the spectra. The second step is similar to the first step: using the new intrinsic 
X-ray spectrum for a new SED, which is in turn used to generate new set of photoionization models. 

The Chandra mid-state data (ID 4399) was observed simultaneously in the UV band by FUSE . We adopted this UV data point 



1 Dunn et al.ll2008f) . while other data were obtained from NASA/IPAC Extragalactic Database (NED): radio data from VLA, 25 fi data 



from IRAS, and r band data from SDSS. 

For each state, the flux point at 0.5 keV band on the intrinsic X-ray continuum was simply linked to the UV point and to other data in 
the SED by straight lines in log-log space, as shown in Fig[5] We derived luminosities of 0.5-13.6 keV band from the intrinsic X-ray spectral 
models, measured luminosities of 13.6 eV - 0.5 keV band from the two point lines in SEDs, and added them to determine Li 0n . Detailed 
intrinsic continuum parameters and luminosities are listed in Table|7] 



1 http://heasarc.gsfc.nasa.gov/docs/software/xstar/xstar.html 
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4.2 Warm Absorbers from Model Results 

4.2.1 One WA Component 

Turner et alj 1 19961) modeled the ASCA spectrum of Mrk 290 with one WA component generated using XSTAR, obtaining an ionization 



parameter log{^) — 1.38 and a column density Nh of 8 x 10 21 cm~ 2 . Similarly, we began with one WA in fitting Chandra HETGS spectra. 
The Chandra spectra were binned to 0.01 A, and fitted with one ionization component. The model represents most absorption lines well, and 
the best fit results are shown in Table[6] The absorber is blue-shifted by —450 ± 30 km s" 1 relative to the systemic redshift, and neither its 
column density nor ionization parameter changed apparently from mid-state to high-state during half a month. Since the UTA feature is quite 
clear in the RGS spectrum, one WA component fitting is not suitable for the low-state data. 



4.2.2 Two Phase Warm Absorber 

Coarsely binned high- and mid-state spectra are shown in Fig. [6] in which the green lines represent the best fitting model from Table[6] The 
visual inspection of the spectra reveals the unfitted UTA feature clearly in the 16-17.5 A band, suggesting that the single WA component is 
not sufficient to account for all absorbing features. 

The UTA is a broad feature, and its shape is sensitive to the charge states of Fe producing it, so we rebinned the spectra heavily. Since the 
UTA is a key diagnostic of material in a low-ionization state, a second WA with a lower ionization level is required to be included in the model. 
We modeled the high-state spectrum using a two phase WA with parameters set free, except the redshift of the higher ionization component 
was fixed to z = 0.0289 from Table|6] The two phase warm absorber provides a good fit to the high-state spectrum as shown in Fig.|6]by the 
red line model. The Fe UTA is well fitted by a lower ionized component with log(^) = 1.62 ± 0.15, Nh = 5.4(±1.8) x 10 20 cm" 2 . The 
outflow velocity is similar to that of the high ionization component, and the other best fit parameters are reported in Table[7] 

Providing that no properties in Table [7] varied in 17 days and the dynamical timescale is much longer than the radiative recombination 
timescale, when fitting the mid-state spectrum, we fixed the redshift and column density of each absorbing component to that of high-state, 
but allowed the ionization parameter to vary. Unfortunately, the large uncertainties in £ of the low ionization component cannot constrain the 
possible response of opacity to the continuum variation. 

The two phase WA model was also used to fit the RGS spectrum, adding several gaussian lines to account for the local absorption lines 
and the O VII forbidden emission line. Fig.[7]shows the best model for RGS spectrum, and fitting results are listed in Table|7] As expected, the 
deep UTA feature comes from the lower ionized material, while other features identified as transitions of Ne IX, Ne X, Fe XVIII, Fe XX and 
N VII, come from the higher ionization component. The C VI absorption lines are believed to be associated with the low ionization component 
), though the observed lines are much deeper than what the photoionization model predicts here. However, it could also be 
the result from carbon-rich gas. The Fe XVIII complex is not well fitted here, as in HETGS spectra (as shown in Fig.[6}, which is partly due to 

i, and partly due to the uncertainty of the iron abundance 
(solar abundance is used in the photoionization models). Neither the ionization parameter nor the column density of the two WA components 
varied significantly compared to that found in the HETGS spectra, but the outflow velocity increased by a factor of about 2.5. 



i Sako et alJbOQl 



the uncertainties of di-electronic recombination rates for iron L-shell ions Jgu]|2003| 



4.2.3 The Third, Hottest Ionization Component 

We have applied the two phase warm absorber model back to the 0.01 A binned high-state spectrum again and checked residuals carefully. 
Three important Fe L-shell absorption lines were left unfitted, Fe XXI (12.284 A), Fe XXII (11.770 A), Fe XXIII (10.981 A). A hotter high 
ionization component is required to explain these absorption lines, as labeled in Fi g.H] Such a very high i onization absor bing component 



has b een reported in a few sources, e.g. NGC 3783 iNetzer et al]|2003h . NGC 5548 Jsteenbrugge et al.ll2005b and NGC 985 jKrongold etal] 
20091) . 



We added the third ionized absorber in the high resolution high-state spectrum to test for the possible presence of this component. The 
lowest ionization component was fixed to the values of the two phase WA model fit, where they were best determined. The fit parameters are 
given in Table[8] The Cash statistic gives a difference of 45 for adding 3 d.o.f (2590/2341) compared to the model with one WA. The model 
of three ionization components are shown in Fig. [8] The comparison of the best-fitting model to the high-state spectrum can be seen in Fig. [5] 
in which the flux spectrum is binned to 0.025 A for clarify. 



4.3 The L-shell silicon lines 

In Chandra HETGS s pectra, Si L-shell absor ption lines do not exist in the fitted model, because so far they have mostly been omitted from 
the XSTAR database ( McKernan et alj|2007 ). Fig. 1101 shows these unfitted lines in the mid- and high-state spectra. The mid-state spectrum 
shows Si XI, Si IX, and one significant line (at 7.23 A) that could be Si VII with an outflow velocity 2000 km s" 1 , Si VIII with an inflow 
velocity 700 km s" 1 , or spurious. In the high-state spectrum Si XI, Si X, and Si IX are evident. 

The 7.23 A line disappeared totally while Si IX got stronge r from the mid -state to the high-state. Furthermore the ratio of Si IX to 
Si XI changed from 0.6 ± 0.5 to 1.6 ± 0.4 over the two states. Netzer 1 20041) sugge sted that the Si L-shell absorption lines should be 
generated from the same ionization component producing the Fe UTA. iGu et al. I J2006h calculated the new Fe M-shell data and confirmed 
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this suggestion, especially for the NGC 3783. The variation of Si L-shell lines in Mrk 290 demonstrated the change of the ionization parameter 
in the lowest ionization component. 



5 DISCUSSION 



5.1 The Location of The WAs 

When the ionizing luminosity Li OTl and the ionization parameter £ are known, one would expect to constrain the distance R by estimating 
the electron density n e , since £ = Li 0n /(n e R 2 ) . Studying the variation based on non-equilibrium models is one method to constrain n e 



l Nicastro et al J 19991) . An excellent approximation for the time to reach a new photoionization equilibrium is 



Otrec , -^e) eq n*e 



[arecix*- 1 ,T e )/a rec (x i ,T e )] eq + [n x ,+i/n x i] 



(1) 



where Xi is the relative density of ion i of the element X, and the recombination time a rec (x' i ,T e ) can be obtained from Shull & van 
Steenberg (1982) data. 

During Chandra observations, the luminosity increased by a factor of 1.4, thus the log(£) was expected to increase 0.15 dex. However, 
the ionization parameter of the medium ionization component, which is best constrained in the Chandra spectra, did not change (Table |6). 
Providing the flux increase happened between June 29 and July 15 in 2003 and the high-state observation lasted for 2.5 days, the timescale 
for the gas to reach the new photoionization equilibrium must be longer than 2.5 days. Considering the ions O VI - O VIII and Ne VIII - Ne X, 
we find the electron density of the medium ionization component n e < 6.2 x 10 4 cm -3 . Combined with Li on and £, the distance from the 
center black hole to this WA is further than 0.9 pc. 

The large error bars of £ in the lowest ionization component make it hard to check whether £ varied. However, based on the variation 
of Si L-shell lines, we conclude that the ionization state of the lowest ionization component changed over the 17 day span. Using the same 
method as above, we find that n e > 5.5 x 10 4 cm -3 and that the distance to this component is less than 2.5 pc. It might be common tha t 
a lower ionization component reaches a new photoionization equilibrium easier than a higher ionization components dKrongold et al ]l2005h . 
as in NGC 3783, NGC 4051, NGC 5548. 



5.2 



On The Origin of The Emission Lines 

The location of gas produ cing soft X-ray emission lines is unclear. These lin es may originate from just the WA l Netzer et"al1 2003h. very 
close to or within the BLR i Costantini et alJ2o"o7 ; Longinotti et alj2008 , 201fj|) . within the Narrow Line Region (NLR) l Guainazzi & Bianchi 



2007), or between the BLR and the NLR l lNucita et al .Il2010f) . The soft X-ray emission lines in Mrk 290 are redshifted with respect to the 



systemic velocity, and the velocity did not change over 3 years, so they ar e probably not related to the W As. 



Vestergaard & PetersorfeOQfjl) , indicating a possible origin in the 



The FWHM of Fe k a is similar to the width of H/3 (5323.5 km s" 
BLR i Bianchi et al.112008 ). However, forcing the lines to be narrow d oes not make the fits significantly worse, which suggests that the torus 
is still the primary origin for this fluorescent line dBianchi et al ]l2007h . 



5.3 The Torus Origin 

With the reverberation mapping results i Dennev et al. 2O10h . the black hole mass of Mrk 290 is 2.43 ± 0.37 x 10 7 M Q . The time delay 
of H/3, 8.72l^Q2 days, indicates that the radius of the BLR is 0. 007 pc. The thermal bolometric luminosity calculated from the SED, 
Lfc ; = 2.5 x 10 44 ergs -1 , gives the radius of the torus: l.lpc (e.g. Krolik & Krissll2001 ). Thus, the radius of the torus seems similar to the 
distance of the WAs. 

The thermal photoionization equilibrium curve i Krolik et al.ll 198 1 > is presented in Fig.QT]for the high-state SED. The curve marks the 
points of thermal equilibria in T - S plane, where T is the electron temperature and the "pressure" ionization parameter H = £/47rcfcT. The 
three WAs best modeled in the high-state spectrum are roughly in pressure balance as shown in Fig.QT| indicating that all three components 
may occupy the sa me spatial region. Com bined with the distance, we suggest that the origin of the WAs is a thermal wind from the torus, the 
model proposed by Krolik & Krissl 1 2001 ). 

The WAs in the RGS spectrum do not show significantly different ionization parameters and column densities compared to that found 
from the HETGS spectra. This supports the model that the WAs come from the torus. The large change of the outflow velocity, from 
450 kms -1 to 1260 kms" 1 , gives a large acceleration 0.6 cms . Radiation pressure is an attractive means for accelerating the gas 



l Crenshaw et al]|2003f) . and the acceleration can be written as a = 



Gil 



where k is the absorption cross-section per unit mass. 



For Mrk 290, the acceleration by radiation pressure is ~ 0.01 cm s~ , much less than the observed one. As a result, an entirely new WA gas 
may have passed through our line of sight during the XM M-Newton observat ion. 

A natural geometry for WA winds is bi- conical I Dorodnits vn et al 1l2008l) . Assuming this, the mass loss rate can be estimated. We use a 
formula for the mass outflow rate derived by Krongold et al. 1 2007h : 
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Mont = 0.87Tm p N H VrRf(S, 0) (2) 

where f(S, <j>) is a factor that depends on the particular orientation of the disk and the wind, and for all reasonable angles (S > 20° and 
4> > 45°) is of order of unity. With the observed v r (the line of sight outflow velocity), Nh, and the distance of torus R = 1.1 pc, the 
mass outflow rate contributed by the medium and lowest ionized WAs is ~ 1-Mq per year, of the same order as the values in the simulation 



l Dorodnitsvn et al. 20081) . The outflow mass rate estimate is larger than the accretion rate, M acc = i^r 1 ~ 0.5Mq per year, assuming 

T] = 0.1. 

The escape velocity (v esc — y/2GM/R) from the location of the torus is about 440 kms -1 , which is nearly identical to the outflow 
velocities of the WAs in Mrk 290. The outflow velocity is only a lower limit to the actual speed, thus the warm gas would not likely fall back, 
and would move 5 kpc away from the central region into the host galaxy over the AGN lifetime (10 7 yr). The total mass lost due to the WAs 



is Mon t ~ 10 7 Mq, and the integrated kinetic energy (given by M ou tv 2 /2) is ~ 2 x 10 5j ergs in 10 7 yr. As suggested bv lKrongold et al 



1 20071) . the hot interstellar medium in the host galaxy can be heated to 10 7 K and would evaporate, ceasing star formation. As a consequence, 



this outflow mass rate may affect the evolution of the host galaxy and the surrounding intergalactic medium, over the AGN lifetime. 



5.4 The Connection between UV and X-ray Absorbers 

Mrk 290 was observed by FUSE simultaneously with Chandra . Kaiser et al. 1 20061) detected three intrinsic absorbers by the O VI doublet 
in the UV band, blueshifted by 526, 470 and 439 km s -1 with respect to the systemic velocity. Comparing them with the outflow velocities 
of the three WAs in X-ray band, 540 ± 150, 450 ± 30 and 390 ± 60 kms -1 , we find the UV absorption components correspond to the 
X-ray WAs within the uncertainties. One geometry to interpret this connection is that high-d ensity, low- i onization UV-absorb ing clouds are 
embedded in a low-density, high-ionization material that dominates the X-ray warm absorber i Krissl20o"4 ). Kaiser et al. 1 2006b estimated that 
the hydrogen column density is 9 x 10 14 cm -2 . Assuming n e > 6 x 10 4 cm -3 , the thickness of the UV-absorbing clouds must be less than 
1.5xl0 5 km. In fact, many X-ray observed absorption comp onents appear to be accompanied by corresponding UV absorption components 
with the same outflow velocity (e.g. Steenbrugge et alj|2009 



6 CONCLUSIONS 

We summarize the main results of this analysis of the Mrk 290 high resolution X-ray spectra as follows: 

1. The Chandra -HETGS spectra reveal complex absorption features including the Fe M-shell unresolved transition array (UTA), indicat- 
ing the presence of intervening ionized absorbing gas. Photoionization modeling shows that a combination of three ionization absorbers with 
similar kinematics (outflow velocities ~ 450 km s" 1 ) can fit the data and are consistent with the three absorption components in the previous 
far UV band study. The ionized component with highest ionization is required by the presence of several absorption lines of Fe XXI-XXIII, 
which are rarely observed. 

2. The ionizing continuum increased 40% during the Chandra observation, on a timescale of 17 days. Using the lack of response in the 
opacity of the higher ionization component to this fluctuation, we constrain the distance from the ionizing source to > 0.9 pc . Si L-shell line 
ratios show significant changes and are related to the UTA, so they result from same gas with the ionization level. If the variation derived 
from the Si L-shell lines of the lowest ionization component is real, we may set an upper limit of 2.5 pc for the distance to the central source. 
The three ionization components lie on the stable branch of the S-curve, showing roughly the same gas pressure, which indicates that they 
occupy the same spatial region. We conclude that the warm absorbing gas originates in the thermal wind from the torus. 

3. During the XMM-Newton observation, the ionizing luminosity of Mrk 290 was only half of the high luminosity of the Chandra 
observation. The RGS spectrum is well fitted by a two phase warm absorber, with several additional local absorption lines attributed to the 
high velocity cloud Complex C. Neither the ionization parameter £ nor the column density Nh of the two absorbing components varied 
significantly, compared to the results from Chandra observations. The outflow velocities of both components increased by a factor of 2.5. 
However, the radiation pressure in Mrk 290 is not enough to accelerate the gas to ~ 1260 kms -1 . We suggest that entirely new WA gas 
from the torus passed through our line of sight. 

4. The estimated mass outflow rate of the ionized wind is ~ IMq per year, larger than the nuclear accretion rate ~ O.5M0 per year, 
and the input kinetic energy is ~ 2 x 10 ergs. 

5. In both high- and low-state spectra the Fe Kq emission lines are detected. The O VII and Ne IX forbidden lines are the most prominent 
soft X-ray emission lines, with a mean redshift velocity of ~ 700 kms -1 with respect to the systematic velocity. There seems to be no 
relation between emission lines and warm absorbers. 
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time [ks] time [ks] 

Figure 1. Light curves from Chandra observations of Mrk 290. Each panel corresponds to an observation in Table 1. The upper half of each panel shows the 
hard X-ray (2-10 A) count rate and the lower half shows the softness ratio, defined as the count rate in the 15-25 A band divided by that in the 2-10 A. The 
dashed lines serve as variability references. The count rate and softness ratio are lowest in observation 4399, while the other three observations are consistent 
with each other. 
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Figure 2. Combined MEG and HEG 1 st order spectra compared to the best-fitting power-law plus black body continuum model. Residuals are given in the 
bottom two panels for the high-state (black line) and for the mid-state (green line). The data were binned to 0.01 A. Possible spectral features are marked. In 
the residuals, C denotes for the Cash Statistic, and V AC is a proxy for \ because AC is distributed as \ 2 ■ 
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Energy (keV) 

Figure 3. Comparison of the XMM-Newton pn spectrum to the best-fitting power-law plus black body model in the 0.5-10 keV band. Iron K« emission and 
an absorbing feature from 0.7 to 0.9 keV can be seen clearly in the ratio of the data to the model. The absorption is due to the warm absorber. 
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Figure 4. Comparison of the RGS spectrum to best-fitting power-law plus black body model from pn spectral fit. The Fe UTA is clearly observed in the 16.0 - 
17.5 A range, along with other features that are marked. 
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Figure 6. The photoionization models for coarsely binned HETGS spectra. The upper panel shows the high-state spectrum and the lower panel shows the mid- 
state spectrum. In the residuals, the green lines indicate the model with one WA component, while the red lines represent the model with two WA components. 
The difference between the two models is primarily found at the Fe UTA feature. The residuals around 23 A are partly due to the O VII(f) emission line. 
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Figure 7. The photoionization model for the low-state RGS spectrum. The lines labeled in black are at z = and are due to Galactic gas. The VII(f) 
emission line is significant. 
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Log(C) = 1.62 Log(N H ) = 20.7 
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Figure 8. Three WA components are used to model the high-state spectrum of Mrk 290, shown on a scale where 1 is the incident continuum level. The three 
lines of Fe XXI-XXIII marked in black indicate the highest ionized component, while the broad Fe UTA results from the lowest ionization component (in 
green). 
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Figure 9. The best fit model for high-state spectra in the 1.5-25 A band, which consists of three WA components. Four fitted emission lines are marked in 
purple. The upper part and lower part of each panel show the count spectra and the flux spectra, respectively. 
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Figure 9. Continued. 
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Figure 10. The 6.7-7.5 A regions of mid- and high-state spectra are shown in the left and right panel, respectively. The Si Vll/Si VIII line at 7.23 A 
significantly stronger in the mid-state than in the high state, indicating that the ionization state of the WA responded to the continuum changes. 
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Figure 11. Thermal stability curve for a low-density gas illuminated by the high-state continuum of Mrk 290. The quantity H = i^/A-KckT is the pressure 
ionization parameter. The black sections show the locations of the three warm absorbers, with the corresponding uncertainties on the ionization parameters. 
The similar value of H indicates that the three components have roughly the same gas pressure. 
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Table 1. Chandra HETGS observation log 



ID 


Start Time 


Exposure time 


4399 


2003.06.29 01:38:04 


85.1 ks 


3567 


2003.07.15 01:18:47 


55.1 ks 


4441 


2003.07.15 22:54:28 


60.9 ks 


4442 


2003.07.17 11:18:15 


50.2 ks 



Table 2. XMM-Newton observation log 



ID 


Start Time 


Filtered time 


0400360201 


2006.04.30 16:47:58 


20 ks 


0400360301 


2006.05.02 15:21:00 


7.5 ks 


0400360601 


2006.05.04 14:45:53 


15 ks 


0400360801 


2006.05.06 14:39:11 


14 ks 



Table 3. Continuum parameters 



Date State T kT Flux2-iokeV L,2-i0kcV IDs and notes 

(eV) (10" 11 orgs- 1 cm~ 2 ) ( 10 43 org a ~ 1 ) 



2003.06.29 high- 1.86 ±0.02 81 ± 6 1.71 3.13 3567,4441,4442 
2003.07.15 mid- 1.55 ± 0.02 69 ± 9 1.20 2.19 4399 

2006.04.30 low- 1.58 ±0.01 87 ± 1 0.88 1.61 EPIC pn, RGS 



"The photon index of the power-law component. 'The temperature of the black-body component. 
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4. Absorption lines from the high resolution spectra with Poisson probability greater than 0.995 





Flux 


Outflow Velocity 


EWs 


Poisson 


Ion Name & A rest 


(A) 


(10~ 6 ph cm~ 2 a -1 ) 


(kma- 1 ) 


OnA) 


Probability 


(A) 


High-state spectrum 


4.202 ± 0.006 


3.7 ± 1.8 


790 ± 420 


7 ± 4 


0.9990 


S XV (4.088) 


6.360 ± 0.002 


2.6 ± 0.5 


530 ± 90 


5 ± 1 


0.9995 


Si XIV Lya (6.183) 


6.844 ± 0.003 


2.5 ± 0.5 


270 ± 110 


5 ± 1 


0.9993 


Si XIII(r) (6.648) 


7.133 ± 0.001 


2.3 ± 0.5 


380 ± 60 


5 ± 1 


1.0000 


Si IX (6.923,6.939) 


8.664 ± 0.003 


3.8 ± 1.4 


440 ± 120 


8 ± 3 


1.0000 


Mg XII Lya (8.421) 


9.433 ± 0.002 


3.4 ± 1.0 


480 ± 80 


8 ± 2 


0.9975 


Mg XI (9.169) 


9.979 ± 0.005 


3.9 ± 1.4 


520 ± 170 


10 ± 3 


0.9996 


Ne X Ly 7 (9.708) 


10.931 ± 0.008 


3.8 ± 1.8 


630 ± 230 


8±4 


0.9969 


Fe XIX (10.630) 


12.109 ± 0.005 


4.2 ± 2.0 


470 ± 110 


9±4 


0.9988 


Fe XXII (11.770) 


12.482 ± 0.003 


10.0 ± 2.4 


520 ± 70 


21 ±7 


1.0000 


NeX Lya (12.134) 


12.632 ± 0.007 


6.8 ± 2.8 


620 ± 170 


17 ±7 


1.0000 


Fe XXI (12.284) 


13.205 ± 0.010 


9.5 ±4.6 


450 ± 240 


23 ± 11 


0.9999 


Fe XX (12.824,12.846,12.864) 


13.841 ± 0.010 


21.0 ± 7.1 


330 ± 220 


45 ± 15 


1.0000 


Ne IX (13.447), Fe XIX (13.423,13.375) 


13.900 ± 0.001 


8.8 ±2.0 


400 ± 30 


19 ±4 


1.0000 


Fe XIX (13.497,13.518) 


13.902 ± 0.008 


8.9 ±5.5 


350 ± 180 


21 ± 13 


1.0000 


Fe XXI (13.507) 


14.613 ±0.006 


10.1 ±4.6 


570 ± 130 


23 ± 11 


1.0000 


Fe XVIII (14.208) 


14.950 ±0.015 


5.1 ±4.2 


530 ± 300 


12 ± 10 


0.9998 


Fe XVIII (14.534) 


15.446 ± 0.006 


7.4 ±4.7 


490 ± 120 


16 ± 10 


1.0000 


Fe XVII (15.014) 


15.621 ± 0.015 


12.1 ±6.3 


330 ± 300 


24 ± 12 


0.9999 


O VIII Ly-y (15.176) 


16.472 ± 0.017 


8.5 ±8.2 


380 ± 320 


17 ± 16 


0.9950 


vm Ly /3 (16.006) 


16.841 ± 0.006 


9.2 ±3.6 


310 ± 110 


18 ±8 


1.0000 


Fex (16.360) 


18.942 ± 0.006 


14.5 ±4.9 


-400 ± 90 


24 ±8 


0.9990 


*0 vm Lya (18.967) 


19.164 ±0.002 


23.4 ±6.2 


460 ± 40 


35 ±9 


0.9994 


VII (18.627) 


19.519 ±0.010 


18.4 ±9.3 


390 ± 170 


27 ± 14 


0.9950 


VIII Lya (18.967) 


Mid-state spectrum 


6.357 ± 0.005 


2.1 ± 1.0 


660 ± 220 


8±4 


0.9983 


Si XIV Lya (6.183) 


7.229 ± 0.002 


2.6 ±0.7 




12 ±3 


1.0000 


Si vm (6.9990) or Si VII (7.063) 


8.665 ± 0.005 


2.3 ±0.7 


430 ±170 


11 ± 3 


1.0000 


Mg XII Lya (8.421) 


10.940 ± 0.005 


2.6 ± 1.3 


370 ± 140 


13 ±6 


0.9987 


Fe XIX (10.630) 


11.348 ±0.006 


2.7 ± 1.2 


440 ± 170 


14 ±6 


0.9978 


Fe XIX (11.029) 


12.481 ± 0.003 


8.6 ±2.5 


540 ± 80 


21 ±6 


0.9958 


Nex Lya (12.134) 


13.205 ± 0.010 


9.5 ±4.5 


450 ± 250 


23 ± 11 


0.9986 


Fe XX (12.824,12.846,12.864) 


15.624 ±0.003 


8.7±3.1 


260 ± 50 


42 ± 15 


0.9959 


vm Ly7 (15.176) 


Low-state spectrum 


12.396 ±0.081 


8.7 ±8.0 


2640 ± 2000 


62 ± 57 


0.9984 


Nex Lya (12.134) 


13.823 ±0.023 


11.8 ±3.3 


1400 ± 520 


73 ± 20 


1.0000 


Ne IX (13.447), Fe XIX (13.497,13.518) 


14.550 ±0.018 


14.0 ±3.1 


1890 ± 380 


83 ± 18 


1.0000 


Fe xvm (14.208) 


21.568 ±0.056 


22.1 ± 11.8 


-490 ± 780 


100 ± 41 


0.9988 


*0 vn(r) (21.602) 


24.833 ± 0.023 


13.8 ±6.6 


-100 ±280 


52 ± 25 


0.9981 


*N VI (24.898), *N VII Lya (24.785) 


25.492 ± 0.029 


11.6 ±6.6 


520 ± 350 


44 ±25 


0.9977 


N VII (24.779,24.785) 



* Galactic absorption lines, the velocities of which are relative to the heliocentric system. 
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Table 5. Emission lines in the high- and low-state spectra 



Ion 


^rest 
(A) 


A obs flux 

(A) (10~ 6 phcm" 2 a -1 ) 


Outflow Velocity 

(kma- 1 ) 


High-state spectrum 


Fe Ka 
Fe xxv 
Ne ix(f) 

O VII(f) 


1.938 
1.859 
13.669 
22.098 


1.995 ±0.008 4.3 ±3.3 
1.897t° Q2g 10.9 ±7.9 
14.115 ±0.007 8.8 ±4.6 
22.831 ±0.008 54.7 ±37.3 


-670 ± 150 
-830 ± 110 


Low-state spectrum 


Fe Ka 
O VII(f) 


1.938 
22.098 


2.001 ±0.021 4.8 ±1.7 
22.818 ±0.026 29.0 ±16.5 


-660 ± 350 



Table 6. Best fit parameters for one WA component 



State 


Log? 


N H 


"Velocity 


Cstat/d.o.f. 


b AC 




(erg, am- 1 ) 


(10 21 cm -2 ) 


(km a" 1 ) 






high- 


2.43 ±0.02 


4.9 ± 0.5 


450 ± 30 


2635/2344 


431 


mid- 


2.41 ±0.07 


4.8 ± 1.3 


450 ± 60 


2071/1843 


92 



a The outflow velocity of WA is relative to the systematic velocity of Mrk 290. 
b Change in the C statistics between models with and without a warm absorber. 



Table 7. The two phase WA parameters for three states. 



State 


°r 


b kT 

(=V) 


Ll-1000Ryd 
(10 43 orgs" 1 ) 




Log? 

(«,.«-'] 


N H 

(10 21 cm -2 ) 


c Velocity 

(kma- 1 ) 


Cstat/d.o.f. 


high- 


1.91 ± 0.01 


119 ± 3 


13.7 


WAl 
WA2 


1.62 ±0.15 
2.45 ± 0.04 


0.54 ± 0.18 
5.33 ± 1.17 


540 ± 150 
d 450 


414/306 


e mid- 


1.64 ± 0.02 




10.0 


WAl 
WA2 


1.52 ±0.55 
2.37 ± 0.08 


0.54 
5.33 


540 
450 


172/152 


low- 


1.69 ±0.01 


105 ± 2 


7.2 


WAl 
WA2 


1.57 ± 0.15 

z - ou -0.08 


1.12 ± 0.18 
3 16+ 035 


1260 ± 160 
1260 ± 250 


225/183 



"The photon index of the power-law component in the intrinsic spectra. 
b The temperature of the black-body component in the intrinsic spectra. 
c The outflow velocity of WA is relative to the systematic velocity of Mrk 290. 
d The outflow velocity is fixed to the value in Table|6] 
e The column density and the outflow velocity of WAs are fixed to the values obtained from the high-state spectral fits. 



Table 8. Best fit parameters for the three WAs in the high-state spectrum 



Log? N H "Velocity 

(crgscm- 1 ) (10 21 cm~ 2 ) (kma -1 ) 



6 WA1 1.62 0.54 540 

WA2 2.42 ±0.04 4.06 ± 0.62 450 ± 30 
WA3 3.20 ±0.14 3.49 ± 1.47 390 ± 60 



a The outflow velocity of WA is relative to the systematic velocity of Mrk 290. 
b The parameters of this lowest ionization component is fixed to the values in Table|7] 



Lagm - 3.20 Log(K H ) = 21.5 
Log(0-2,d-2 Log(Nu}-21.6 

LogfO - 1.62 Log(Nn) - 2C.7 



i 




